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Confexf. Maser transitions are commonly observed in media exhibiting a large range of densities and temperatures. They can be used 
to obtain information on the dynamics and physical conditions of the observed regions. In order to obtain reliable constraints on the 
physical conditions prevailing in the masing regions, it is necessary to model the excitation mechanisms of the energy levels of the 
observed molecules. 

Aims. We present a numerical method that enables us to obtain self-consistent solutions for both the statistical equilibrium and 
radiative transfer equations. 

Methods. Using the standard maser theory, the method of Short Characteristics is extended to obtain the solution of the integro- 
differential radiative transfer equation, appropriate to the case of intense masing lines. 

Results. We have applied our method to the maser lines of the H 2 molecule and we compare with the results obtained with a less 
accurate approach. In the regime of large maser opacities we find large differences in the intensity of the maser lines that could be as 
high as several orders of magnitude. The comparison between the two methods shows, however, that the effect on the thermal lines 
is modest. Finally, the effect introduced by rate coefficients on the prediction of H2O masing lines and opacities is discussed, making 
use of various sets of rate coefficients involving He, o-H 2 and p-H 2 . We find that the masing nature of a line is not affected by the 
selected collisional rates. However, from one set to the other the modelled line opacities and intensities can vaiy by up to a factor ~ 2 
and ~10 respectively. 



1. Introduction 



The first maser signal was detected in 1965 dWeaver et al.l 
119651) . It was due to OH and the intriguing nature of this de- 
tection was r apidly attributed to a non-thermal phenomenon 
(Davies et al. L [196% which originates in the amplification of 
the radiation due to stimulated emission. Since then, many 
improvements have been made both in observational facili- 
ties and maser theory, leading to an abundant literature on 
the subject. The main results, both t heoretical and ob serva- 
tional, have been summarised by, e.g., lElitzurf (1 19921) or iGravl 
(1999). From the observational point of view, major improve- 
ments have been obtained in the last 20 years thanks to the 
rise of sensitive interferomete rs and VLBI techniques (see e.g. 
Humphreys, 2007; Fish, 2007, and references therein). The sizes 
of the regions where maser emission is detected have typi- 
cal scales of a few AU. Nevertheless, some masers arise from 
large spatial scales, such as th e 183.3 GHz line of water vapour 
dCernicharo et all Il990i Il994l) . From the point of view of the 
theory, despite the fact that many improvements have been made 
in explaining the pumping mechanisms of molecular masers, 
there are no definitive explanations that could account for all 
the observations as yet. As an example, the problem of water 
megamasers ob served in the direct i on of AGN is still a sub- 
ject of debate dElitzur eTall 119891: iDeeuchi & Watsonl Il989t 
iKvlafis & Pavlakisll 1999tlCernicharo et al.Ll2006al) . 

The basics of the standard maser theory which is wide ly used 
in ord er to interpret masing lines have been reviewed by lElitzud 
(119921) . The main aspect of this theory is to introduce the effect of 
saturation which puts limits on the amount of energy that can be 
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radiat ed by a masing line. As outlined bv ICollison & Nedoluhal 
(fl995l) . maser theoretical studies adopt two different strategies: 

i) they solve exactly the radiative transfer in the maser but with 
a simplified description of the energy levels of the molecule; 

ii) or deal with an exact description of the molecule but using 
approximations in the radiative transfer treatment. In the latter 
case, the escape probability formalism is widely used, despite 
the fact that the validity of this approximation is questionable in 
the case of masing lines. Indeed, intense masers require a large 
amplification path and, therefore, coherence in velocity over rel- 
atively large scales. In contrast, the basic principle of the LVG 
approximation is that the different regions of the cloud are de- 
coupled due to Doppler shifting, hence limiting the coherence 
path. In this work we present a self-consistent method that aims 
to go beyond the approximations usually made, i.e. we present 
a method that enables us to solve exactly the radiative transfer 
problem and in which we deal with a rigorous description of the 
molecule. 



The theory relevant to the description of masing lines as well 
as the way the radiative transfer is solved is described in Section 
[2] An application for the 0-H2O molecule is presented in Section 
[3]where the method developed in this study is compared to other 
less accurate approaches. Finally, the current method is used in 
Section [4] in order to predict which lines of 0-H2O or P-H2O 
could be inverted for densities and temperatures typical of AGB 
circumstellar envelopes. In this section, the predictions based on 
various sets of collisional rate coefficients involving He, P-H2 or 
0-H2 are also compared. 
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2. Theory 

2. 1 . Absorption coefficient 

In the case of intense masers, the radiation field starts to influ- 
ence the inversion of the population, avoiding a complete fre- 
quency redistribution, leading to differing line profiles for ab- 
sorption and emission. This occurs when the intensity in the 
masing lines exceeds a certain threshold which depends on the 
given physical conditions. In this case, the maser is said to be 
saturated. The way this effect is introduced in the theory is de- 
scribed later in this section. 

Moreover, under a particular geometry the amplification path 
varies with direction, implying that the distribution of molecu- 
lar velocities should in principle b e considered as anisotropic . 
This problem was first identified bv lBettwieser & Kegell(ll974l) . 
In the case of spherical geometry, a common approxima- 
tion, known as "standard theory", consists of neglecting the 
anisotropy of the molecular velocities resulting from the exis- 
tence of directions of different path lengths. The vali dity of this 
assumption was tested in various studies (see e . g. iBettwieserj, 
ll97fllNeufeldll992tlEWering & Watsorlll994l:Eitzuii ll994) 
for spherical geometry. They concluded that considering explic- 
itly the anisotropies of the absorption coefficient could lead to 
variations in the properties of the emerging maser radiation, in 
the case of saturated masers. However, it was pointed out that 
these variations remain unimportant in comparison to the un- 
certainties that affect maser theory, and that the standard theory 
should, at least, be a reasonably good first approximation. Given 
the large number of simplifications adopted by the standard the- 
ory in the treatment of the radiative transfer (RT), we assume 
in the following that the masing line absorption coefficients are 
isotropic, and we follow a p rocedure similar to that described by 
lAnderson & Watsonl (I1993I) . First, we consider the frequency- 
dependent statistical equilibrium equations (SEE) for the level 
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In these equations, Nj represents the total population of level 
j, rij(v) is the population of level j at the velocity v. The quantity 
<p{v) corresponds to the usual Doppler line profile that accounts 
for both the thermal and turbulent motion of the molecules. Ajj 
and Bij are the Einstein coefficients for spontaneous and induced 
emission, Ej is the energy of level j, Cy, are the collisional rates 
from level j to level i and I v (£l) is the specific intensity in a 
direction specified by Q. We note that in equation Q] the elas- 
tic collisional rates Cu have to be formally included since they 
contr ibute to the relaxation of velocities (lAnderson & Watsonl 
Il993l) . Equation[T]can be re-written so that all levels except those 
involving the maser transition are included in the so called pump 
and loss rates. This leads to a phenomenological treatment where 
the full SEE is artificially reduced to a two level system. In the 
following, the labels u and I indicate the upper and lower maser 
levels, respectively. Assuming that for all the other transitions 



the absorption and emission line profiles are identical and are 
given by the Doppler line profile <p{v), we obtain from eq. [T] 

= A„(v) 0(v) - n u (v) [F„(v) + C ul + C uu ] 
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where the pump and loss rates, noted A, and F, (with i e {u, I}), 
are defined by 
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The source function of the masing line is given by: 
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where n„(v) can be obtained from eqs[3]and[7] Contrary to the 
case of thermal lines, the source function is frequency dependent 
for the maser line. 

In order to obtain the source functions and the absorption co- 
efficients from these expressions, we see that the relevant quan- 
tities are the specific line intensities (averaged over the angles), 
and the total level populations, Nj. The latter are obtained by 
solving the frequency independent SEE expressed by: 



= Z [ n j c ;< - NiQ j] + Z [ n j J * - Ni M 
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In the latter expression, the velocity profile (pk(v) = ra*(v)/iV* 
accounts for the departure from the Doppler line profile if the 
line k — » I is a masing transition, for consistency with eq. Q] 
In the case of thermal lines, these profiles are simply given by 
the Doppler profile and are identical for both levels. In this lat- 
ter case, the average of the specific intensities over angles and 
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velocities in eq. [13] can be factorized, which is crucial in es- 
tablishing the precondition ed form of the SEE, as discussed by 
iRvbicki & Hummeil d!991l) . In the case of masing lines and for 
consistency with eq.Q] we have introduced unequal velocity pro- 
files for the maser levels; consequently the factorisation leading 
to the preconditioned SEE is no longer feasible. Hence, the usual 
ALI method cannot be applied. 

In most studies treating the radiative transfer in the presence 
of masers, the pumping rates are taken to be independent of the 
velocity. In this case the expression for the population difference 
is similar to eq. Q except that the velocity dependence of «o(v) 
and J{v) is omitted. Hence, the absorption coefficient will just 
depend on velocity through the term Ri u (y). In the case of un- 
saturated masers, this term is always negligible in comparison 
to J(v) and the absorption coefficient is given by the Doppler 
line profile. However, for a saturated maser, i.e., for a line with 
Riu(v) ~ J> the absorption coefficient decreases at the line centre. 
This translates to the fact that, for a given pumping scheme, the 
molecules pumped to the upper level are readily depleted to the 
lower by stimulated emission since the number of masing pho- 
tons is large. At this stage, the maser amplification starts to be 
linear at line centre and the maser line re-broadens. 



2.2. Solution of the transfer equation 

To solve the transfer equation we have used the short cha racter- 
istic method (SC) introduced bv lOlson & Kunaszl d 1987b . The 
numerical code based on this method has previously been de- 
scribed in iDaniel & Cernicharol ((2008) for the case of thermal 
transitions in ID spherical geometry. Briefly, in the SC method 
we assume that the source function can be reproduced numeri- 
cally by a polynomial function of the opacity in the line (usually 
of order 1 or 2). At each iteration the absorption coefficients are 
known from the current populations, and the opacities are eval- 
uated assuming that they behave linearly between consecutive 
grid points. These assumptions permit us to solve analytically 
the transfer equation at each spatial grid point. In our models, 
the radiative transfer deals with both maser and thermal lines. 
However, the treatment for maser propagation is only switched 
on when the opacity of a line (i.e. t„/ = Ni - gi/g u N„), at a par- 
ticular grid point, is negative. For the same line and for the radii 
where the opacity is positive, the line is treated in a standard way 
relevant to thermal radiation. 

In what follows, we describe a method that allows us to solve 
the transfer equation for masing lines. It is based on the assump- 
tion that the source functions are linearly interpolated. Hence, 
along a characteristic, specific intensities are derived according 
to: 



Ii = 11-10-"*+ an Si-! +/3iSi 



(15) 



where the a, and fij coefficients depend on the opacity A r,- be- 
tween points i and i — 1 . 

In the case of masing lines, a difficulty arises from the fact 
that the absorption coefficient at point i depends on the angular 
average Ru,(y) of the specific intensities /,(#): 
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where the quantities indexed by a cross (e.g. / + ) are obtained 
with the current estimate for the populations of the energy levels 
(i.e. the populations obtained at the previous iteration). The same 
applies to the source functions (see eq. [T2l . Therefore, the ra- 
diative transfer equation becomes integro-differential, whereas 



in the case of thermal lines, the problem consists in solving a 
differential equation. Note that the integral part of this integro- 
differential equation involves an integration of the intensities 
over An steradians. In order to obtain intensities which are fully 
self-consistent with the level populations used at a given itera- 
tion, we thus proceed in three steps. 

2.2.1. First iterative scheme 

As a first step, we derive specific intensities at each grid point 
according to eq. Q3] assuming that the source functions are in- 
dependent of the radiation field. In other words, we use eq. Q~2] 
in order to express the source function assuming that the aver- 
age radiation field is Rt(v). For each frequency, we obtain the 
solution iteratively by calculating I v (6) for 6 e [0 : n] and then 
updating k/ u (v) at grid point i. We iterate until the calculated spe- 
cific intensities are consistent with the initial averaged radiation 
field. 

In the following discussion, Ri„(v) is indexed by i or / when- 
ever it corresponds to the radiation field used in the evaluation 
of Ki u {v) or to the one which results from the evaluation of eq. [15] 
and made use of R' lu (v). Additionally, the solutions obtained dur- 
ing the previous steps are indexed with the superscript > if they 
satisfy R{ u {v) > R\ (v) and with the superscript < on the contrary. 
Once the iterative process has been initialised, a trial value for 
Ri u (v) is obtained according to: 



with 
a = log 



1 - a 



b = log*£(v)-alog«£(v) 



(17) 
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The radiation field derived in this way is used to update ac; m (v). 
Depending on the resulting specific intensities, one of the pair 
(Rf, Ri) or (Rfa, Rfa) is updated. Finally, the iterative process 
is stopped when we have R' lu = R^ , within a given accuracy. 

Figure Q] shows the typical behaviour of r{ as a function of R' lu . 
The averaged radiation fields determined using our method are 
marked with squares and are indexed according to the iteration 
number. 



2.2.2. Second iterative scheme 

Once this process has been achieved, we perform an update of 
the source function at point i according to eq. [12] and using 
the average radiation field Ri u (v) derived during the first itera- 
tive scheme. The whole process is then repeated until we obtain 
convergence for the source function at the current grid point. 



2.2.3. Third iterative scheme 

As previously stated, in order to derive the absorption coeffi- 
cient, it is necessary to know the angular average of specific in- 
tensities over An sr. Using the SC method, the spatial grid is first 
swept from the outermost grid point (i = N) to the innermost 
one (i = 1) providing the intensities /,(#), with i e [1;N] and 
8 e [n/2; n]). The grid is then swept in the other direction so that 
the intensities for 6 e [0;n/2]) are then known. This sweeping 
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10" 10 ' 10 
initial R ]u (i / ) 

Fig. 1. Variation of the averaged radiation field (solid line) as a 
function of the initial radiation field entering in the determina- 
tion of the absorption coefficient. The dashed line corresponds to 
a straight line of slope unity. The solution of the problem corre- 
sponds to the intersection between the two curves. The squares 
correspond to the averaged radiation field obtained at each itera- 
tion. 



unknown inco 
intensifies 




inward propagation 



Fig. 2. Schematic diagram showing the availability of incoming 
intensities during the inward propagation. 



procedure is adopted so that at every moment during the prop- 
agation, only the populations used at the current iteration are 
needed to derive the intensities. In the case of masing lines, an 
additional complexity arises from the fact that during the inward 
propagation, the incoming intensities for 9 < Q\- m have not yet 
been calculated (c.f. figure[2]i with the current populations while 
they enter in the evaluation of the absorption coefficient. 

The third iterative scheme thus consists of obtaining an over- 
all convergence for the intensities /,•(#), with i e [1;N] and 
6 e [Q;tt]. This is done by performing successive sweeping of 
the spatial grid, each iteration leading to an update of the inten- 
sities /;(#). 



3. Application 

In order to test our method, we have performed model calcu- 
lations for the ortho state of the water molecule and the com- 
parison is made with respect to a simplified treatment that aims 
to describe the masing lines. This simplified treatment relies on 
the assumption that the calculation of maser inversion and the re- 
sulting intensities can be achieved in a two-step procedure. First, 
the populations of the molecular levels are determined neglect- 
ing the specificity of maser radiative transfer. Secondly, the in- 
tensities of the masing lines are calculated using the populations 



obtained during the first step and taking into account maser prop- 
agation. This procedure has been use d in previous studies (e.g. 
lYates et allll997l : IWatson et al.Ll2002f) and the way it is presently 
implemented is described in the next sections. 

The models consist of clouds of uniform H2 density, tem- 
perature (assumed to be 400 K in all models), and 0-H2O abun- 
dance (40 levels have been considered). We have adopted a di- 
ameter for the cloud of 4.10 14 cm (2.6 10~ 4 pc) which is the 
typical size of bright maser spots in the 22 GHz line of water. 
Moreover, this size provides a maxi mum amplification path sim- 
ilar to that of the slab depths used in lYates et al.l(11997t) . The col- 
lisional rate coeffic ients are for the collisional system H2O / He 
dGreen et al.L [T993) and are corrected in order to account for the 
different reduced masses of the H2O / H2 and H2 O / He systems. 
We ha ve s elected these rates r ather than those of Dub ernet et al.l 
(200 91) or [Faure et al.l (120071) as they were the ones used by 
I Yates et all d 1997b . These parameters being fixed, we performed 
the calculations in the (n(H2),X(H20)) plane of parameter space, 
with respective values in the range [210 5 ;8 10 9 ] cm -3 and 
[2 10~ 6 ; 8 10 ]. The accuracy of the results will depend on the 
gridding adopted. Typically, convergence is reached when the 
spatial sampling is such that the opacity between two consec- 
utive grid points is of the order or less than 1. For all the cal- 
culations presented here, we adopted a grid with 500 spatial 
grid points, which might be inadequate for the highest densi- 
ties/abundances considered here. However, since we are inter- 
ested in masing lines which are collisionally quenched at these 
densities/abundances, the effect of limited convergence accuracy 
does not affect our conclusions. Moreover, the main goal is to 
compare the effect of the radiative transfer and both model cal- 
culations are performed with the same spatial grid. 

As presented in Section|2] elastic collisional rate coefficients 
for the levels involved in masing lines enter in the definition of 
the absorption coefficients a nd source functions. Test ca lcula- 
tions have been performed by Anderson & Watson] (1 19931) in or- 
der to test the influence of these rates and they concluded that in 
most cases their inclusion produces rather small effects. Hence, 
in the present test calculations, we have assumed that the elastic 
rate coefficients are set to zero. In other words, in the following 
calculations we assume that C„ = in Eq. ?? 

3. 1 . Pumping scheme with/without maser propagation 

As a first step, we compare the results of the exact calculation to 
a simplified treatment that consists of neglecting the amplifica- 
tion inherent to maser propagation when deriving the population 
of the energy levels of the molecule. To do so, we take the ab- 
solute value of the level population difference which then enters 
in the definition of the absorption coefficients and source func- 
tions. These quantities are defined with the usual expressions 
which are suitable for thermal lines. This permit us to treat mas- 
ing lines in the same manner as thermal lines. This approach was 
previously adopte d in plane parallel non-local calculations by 
I Yates et all d 1997b with the aim of predicting which water lines 
could show maser excitation. Their results a re discussed in te rm 
of the gain in the masing lines (see eq. 3 of lYates et all 1 19971) 



7«/(v) = 



A»(v)A M /c 2 
8ttv 2 



(19) 



which is subsequently averaged over the slab depths (the so- 
called depth-wei ghted gain coefficient defined by eq. 7 of 
I Yates et all 1 19971) . Since y„i{v) is proportional to the absorption 
coefficient ki„(v), it is equivalent to consider either the depth- 
weighted gain coefficient or the opacity in the masing line, and 
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in the following, we choose to discuss the results in terms of 
line opacities. We emphasise the fact that at this stage, the ap- 
proximate method described here does not involve any treatment 
specific to maser radiation and that could account for saturation 
effects. Thus, the dependence on frequency of the population dif- 
ference An(v) which enters in the definition of y u i(v) [or A7„(v)] 
is just given by the Doppler line profile <p(y). 



off on 




log[n(H 2 )] log[n(H 3 )] 



Fig. 3. Isocontours of the 0-H2O opacities at line centre for the 
6 16 - 5 2 3 line at 22 GHz, the 4 14 - 3 2 i line at 380 GHz and the 
4 2 3 - line at 448 GHz lines. The left column corresponds to 
the models where the treatment of maser propagation is omitted 
and the right one to the current results (see text). The isocon- 
tour values as well as the correspondence between grey scale and 
opacity values are displayed on top of the two columns. The iso- 
contours that correspond to opacities of -4 and -40 are displayed 
in thick blue lines. Negative values for the opacity correspond to 
red lines. 

Figure |3]shows the opacity at line centre of the 22 GHz (6 16- 
5 23 ), 380 GHz (4i4-3 2 i) and 448 GHz (4 2 3-3 30 ) masing lines, 
obtained in the two treatments (labelled respectively on or off 
according to whether maser propagation is accounted for or not). 
Figure |4] shows the line opacities for the same lines with respect 
to cuts in density or H 2 abundance. We can see on these figures 
that including the current maser propagation treatment leads to a 
substantial decrease in the predicted line opacities for these lines. 
Indeed, without treating the exact radiation field in the masing 
lines, we obtain minimum optical depths of the order of -70 for 
the 22 GHz line and -40 for the 380 and 448 GHz transitions. 
With our method, the minimum optical depth for the 22 GHz 
line is reduced to -20 and to -5 for the other two lines. 



n(H 2 ) = 2 10 7 cm 3 X(H 2 0) = 4 10 




mil T j -60 b m±m lIj a 

icr 5 icr 4 icr 3 10 6 10' 



X(H 2 0) n(H 2 ) 

Fig. 4. Comparison between the opacities derived with (solid 
lines) or without (dashed lines) accounting for maser propaga- 
tion for the masers at 22 GHz, 380 GHz and 448 GHz. The left 
column corresponds to a cut in density at n(H 2 ) = 2 10 7 crrT 3 and 
the differences are shown with respect to variations in X(H 2 0). 
The right column corresponds to a fixed water abundance, i.e. 
X(H 2 Q) = 4 IO- 4 . 



22 GHz 

n(H 2 ) = 4 10 7 cm 




10 5 10 4 10 

X(H 2 0) 



Fig. 5. Brightness temperatures for the 22 GHz line obtained 
with the exact treatment (solid lines) and with the approximated 
treatment (dashed lines). The intensities correspond to a line of 
sight that crosses the centre of the sphere and are given for the 
rest frequency of the line. 



3.2. Two-level maser propagation 

In order to infer the influence of these changes in line opaci- 
ties, we now consider the resulting masing line intensities. To 
do so, any reasonable approximated treatment should include a 
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V LSR t km s '] 



Fig. 6. Emerging brightness temperature of the 22 GHz line for 
the model with n(H 2 ) = 4 10 7 and X(H 2 0) = 6 1(T 4 . The solid 
line corresponds to the intensity predicted by the exact calcula- 
tion and the dashed line to the one predicted by the approximated 
RT (see text). 



RT theory that accounts for saturation effects. To obtain a refer- 
ence point for the current results, we proceed as follows: from 
the populations previously obtained in the approximate treat- 
ment, we compute the pump and loss rates A,(v) and F,(v). These 
values enter in the definition of the population difference An(v) 
through eq. Q Then, for a given masing line, the intensities are 
computed by adjusting Ri u (v) as described in section|2] 

Using this approximation, we obtain line intensities that dif- 
fer largely from those obtained with the exact treatment de- 
scribed in section [2] This is illustrated in Fig.[5]for the 22 GHz 
line. This figure shows the brightness temperatures obtained in 
the two treatments at a density n(H2) = 4 10 7 cm 4 , for the centre 
of the line and for a line of sight that passes through the centre of 
the sphere. This figure shows that for water abundances X(H20) 
> 5 10~ 5 the approximate treatment leads to intensities larger by 
two orders of magnitude. This originates from the opacity over- 
estimation obtained in the approximate treatment which has been 
discussed in the previous section. Additionally, Fig.|6]compares 
the line profiles computed in the two treatments for the model 
with n(H 2 ) = 4 10 7 cm 3 and X(H 2 0) = 6 1(T 4 . In this figure, we 
see that the higher intensity found in the approximate treatment 
is accompanied by a broader line profile. This result is produced 
by the higher line saturation obtained in the approximate treat- 
ment with respect to the exact treatment. This fact is illustrated 
in figure [7] where the ratio Ri u (v)/J(v) is plotted as a function 
of the velocity offset to line centre and of the normalised radius 
in the cloud for the 22 GHz line and for the model with n(H 2 ) = 
4 10 7 crrT 3 and X(H 2 0) = 6 1CT 4 , for both the exact and approxi- 
mated methods. We see in this figure that in the exact calculation 
the degree of saturation is low, i.e. Ri u (v)/J(v) < 1 at any radii. 
In contrast, the saturation is high in the approximate calculation 
with Ri u (v)/J(v) > 10 at most radii. 




(a) exact treatment 




(b) approximated treatment 

Fig. 7. Ratio R[ u (v)/J(v) as a function of the offset to line centre 
and of the normalized radius in the cloud. This ratio is given for 
the 22 GHz line and corresponds to the model with n(H 2 ) = 4 
10 7 and X(H 2 0) = 6 1CT 4 . The two plots correspond to (a) the 
exact treatment and (b) an approximated treatment (see text). 

In the exact treatment, since the masing line intensities are 
taken into account when deriving the populations, the population 
difference between the upper and lower level of the masing line 
tends to be reduced by induced radiative de-excitation. On the 
other hand, in the approximate treatment, the averaged radiation 
field is largely underestimated. This entails that the amount of 
induced radiative de-excitation events is underestimated, which 
leads to a greater difference between the populations of the lower 
and upper levels. Consequently, the derived masing line opac- 
ity is larger which is accompanied by a larger degree of satura- 
tion of the line. So, the differences encountered in the two treat- 
ments originate in the overall flow of population from the upper 
to lower level, which is induced by radiative de-excitation. The 
evaluation of this flow of population differs in the two treatments 
since it relies on the accuracy of the determination of the inten- 
sities in the masing line. 

We note that in section [2] eq.Qis derived from the velocity 
dependent SEE and from this expression it is seen that an in- 
crease in the average radiation field leads to a decrease in An(v). 
From similar considerations dealing with the velocity indepen- 
dent SEE, it can be seen that increasing the value of the radiation 
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Fig. 8. The left column corresponds to the opacity of three purely 
thermal lines which are obtained in a treatment where masing 
lines are treated explicitly. The right column corresponds to the 
ratio r T - t°M /t°" of the opacities obtained with the two treat- 
ments (see text). 



field for the maser transitions would lead to lower population 
differences for the levels involved in the maser transitions. This 
behaviour is not a specific case and we can infer from it that 
considering explicitly maser propagation will necessarily lead 
to a decrease in line optical depths in comparison to treatments 
where the average radiation field would be underestimated. 

3.3. Thermal lines 

The two treatments for maser propagation differ in the evaluation 
of the average radiation field of the masing lines. Considering 
the current maser propagation treatment leads to higher values 
for the average radiation field which, in turn, when inverting the 
SEE, leads to a diminution in the population differences between 
levels involved in maser transitions. The modification of the pop- 
ulations of these levels subsequently modifies the whole energy 
level populations. This is illustrated in figure [8] for three lines 
which only show purely thermal emission. In this figure, we see 
that the opacities of these lines are modified and can show vari- 
ations up to ~ 25% between the two treatments in the region of 
the (n(H2),X(H20)) plane where maser transitions are the most 
opaque. Nevertheless, for the range of parameters considered 
here, the number of lines which suffer from optical depth varia- 
tions greater than 10% is low, and for the majority of the thermal 
lines the opacities agree within a few percent between the two 
treatments. 



4. Maser prediction 

Most of the information that can be extracted from water line 
observations relies on the modelling of its excitation. From 
this point of view, water is a difficult molecule to treat since 
its high dipole moment causes most of its transitions to be 
sub-thermally excited (see, e.g. ICernicharo et al.l l2006bl). har- 
bouri ng very large opacities (see, e.g. Gonzalez-Alf onso et al 
1998), and many of them being maser in nature (ICh eung et al 
1969[ IWaters et all 119801; iPhillips et al. L Il980t iMenten et all 



1990a b; Menten & Mel nick L 1 199 ll Tcernicharo et al 

all" 



1994li l 996L ll999. 2006b; Gonzalez- Alfons o et all 1 1 9951 fl 998) 



1990 



Accurate modeling thus requires, in addition to a good descrip- 
tion of the source structure, the availability of collisional rate 
coefficients. Moreover, as discussed above, special formalisms 
have to be developed in order to have a realistic treatment of the 
transfer of radiation in the case of saturated masers. The method 
developed in the present work for masing lines is used to as- 
sess which water lines are expected to show maser excitation as 
well as the relative opacity of the masers. To do so, we have 
performed calculations for various gas kinetic temperatures, H2 
and water volume densities. The grid was set up using the bound- 
aries : 100 K < T k < 1000 K, 10 6 cm" 3 < n(H 2 ) < 10 10 cm 3 and 
n(H 2 0) e {lO 2 ; 10 3 ; 10 4 ; 10 5 } cm 4 . Additionally, the parameter 
space is reduced by considering that the water abundance can- 
not exceed a relative abundance of 10~ 4 with respect to H 2 . This 
range of values for the considered parameters is appropriate to 
the case of the circumstellar envelopes of evolved stars. The ra- 
tios of the line opacities with respect to the maximum 22 GHz 
opacity are given in Table [TJfor o-H 2 (D and Table[2]for p-H 2 0. 
For the two water symmetries, we compare the o pacities ob- 
tained with the rate coefficients of Faur e et al.1 (120071) (inde xed as 
SET1), iDubernet et all (j2006l 120091) ; iDaniel et al.l (I20T0I I20TTI) 
(indexed as SET2) and Green et all (11993 ). In these tables, the 



ratios obtained with the approximation introduced bv lYates et al.l 
(1997) are indicated in parentheses. 

The parameter space con sidered here is a sub-space of the 
parameter space explored by Yate s et al.l (119971) . since the maxi- 
mum gas temperature in their calculation is 2000 K and because 
they include the pumping by dust photons. In the present case 
we use a lower boundary for the temperature since, in order to go 
higher, we would have to extrapolate the rate coefficients. Since 
we aim to compare various sets of collisional rate coefficients, 
we choose not to introduce artefacts due to the extrapolation pro- 
cedure in the comparison of the results. Dust can play an impor- 
tant role in deriving the populations of water levels. However, 
since we report the ratio of the maximum opacity of the masing 
lines with respect to the maximum opacity of the 22 GHz line, 
the consideration of dust can only affect the current result for the 
lines which are radiatively pumped. For these lines, the ratios re- 
ported in Tables Q] and [2] have to be considered as lower limits. 
On the other hand, the ratios reported should be accurate for the 
lines which are collisionally pumped. 

In or der to compare the current results with the results re- 
ported bv lYates et al.l dl997l) . we have to convert the population 
density ratio of the latter study into opacity ratios. To do so, we 
start from the definition of the opacity along a ray of length L 
along which the density populations are uniform. Assuming that 
the velocity profile <p{v) can be expressed as a gaussian, the opac- 
ity at the line i centre is : 



T(/) 



1 

8^1 



rA(i) g(i) dn(i) xLx 



2k B Ti 



(20) 
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Table 1. Prediction of maser lines for different collisional rate coefficients sets. 



Transition 


Frequency 


He 


P-H 2 


p-H 2 


o-H 2 


o-H 2 






(GHz) 


Green 


SET1 


SET2 


SET1 


SET2 


|t 22 | 






8.2 (10.7) 


5.3 (7.1) 


9.1 (15.1) 


4.9 (6.8) 


7.6 (10.1) 


61,6 - 


5 2 ,3 


22 


1.00 (1.00) 


1.00 (1.00) 


1.00 (1.00) 


1.00 (1.00) 


1.00 (1.00) 


10 2 , 9 


- 9 3 , 6 


321 


0.73 (0.88) 


1.11 (1.16) 


0.67 (0.70) 


1.17 (1.15) 


0.78 (0.89) 


4i, 4 - 


3 2 ,i 


380 


0.53 (1.78) 


0.83 (2.14) 


0.50 (1.85) 


0.93 (2.35) 


0.59 (1.63) 


64,3 - 


55,0 


439 


0.48 (0.63) 


0.83 (1.03) 


0.48 (0.69) 


0.91 (1.05) 


0.60 (0.93) 


7 5 ,2 - 


66,1 


443 


0.02 (0.05) 


0.15 (0.21) 


0.11 (0.13) 


0.16 (0.24) 


0.17 (0.25) 


4 2 ,3 - 


3 3 ,o 


448 


0.42 (1.49) 


0.68 (1.89) 


0.39 (1.69) 


0.74 (1.84) 


0.49 (1.55) 


li,o - 


lo.i 


557 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


53, 2 - 


4 4,1 


621 


0.45 (1.69) 


0.76 (2.31) 


0.42 (1.77) 


0.82 (2.13) 


0.53 (1.84) 


7 2 ,s - 


81,8 


1146 


N (N) 


< 0.01 


N (N) 


< 0.01 


N (N) 


3i, 2 - 


2 2 ,i 


1153 


0.02 (0.02) 


0.17 (0.18) 


< 0.01 


0.20 (0.26) 


< 0.01 


63,4 - 


54,1 


1158 


0.18 (0.80) 


0.36 (1.61) 


0.15 (0.56) 


0.44 (2.07) 


0.24 (1.37) 


85,4 - 


7 6 ,i 


1165 


0.09 (0.11) 


0.21 (0.35) 


0.04 (0.04) 


0.21 (0.37) 


0.15 (0.24) 


74,3 - 


6s, 2 


1278 


0.17 (0.49) 


0.38 (1.50) 


0.14 (0.41) 


0.43 (1.70) 


0.21 (0.79) 


8 2 ,7 — 


7 3,4 


1296 


0.31 (0.76) 


0.54 (1.70) 


0.30 (0.90) 


0.59 (1.89) 


0.34 (0.78) 


84,5 - 


9i,s 


1308 


N (N) 


N (N) 


N (N) 


N (N) 


N (N) 


6 2 ,5 - 


33, 2 


1322 


0.19 (0.34) 


0.41 (0.95) 


0.21 (0.27) 


0.46 (0.99) 


0.26 (0.45) 


85,4 - 


%! 


1596 


N (N) 


N (N) 


N (N) 


N (N) 


N (N) 


84,5 - 


7 5>2 


1885 


0.01 (0.02) 


0.06 (0.06) 


0.01 (0.01) 


0.10 (0.10) 


0.03 (0.03) 


5 2 , 3 - 


4 3 , 2 


1918 


< 0.01 


< 0.01 


< 0.01 


0.01 (0.01) 


<0.01 


83,6 - 


74,3 


2245 


< 0.01 


0.20 (0.30) 


0.01 (0.02) 


0.27 (0.43) 


0.01 (0.02) 


73,4- 


64,3 


2567 


< 0.01 


< 0.01 


N(N) 


<0.01 


< 0.01 


83,6 - 




2577 


N(N) 


N(N) 


(N) 


N(N) 


N(N) 


9 2 7 _ 


lOi.m 


2619 


N(N) 


N(N) 


(N) 


N(N) 


N(N) 


9 4 , 5 - 


85,4 


3150 


N(N) 


N(N) 


(N) 


N(N) 


N(N) 



Notes. The first line gives the maximum 22 GHz opacity (t 22 ) obtained for the density and temperature range covered in this study (see Sec.[4]l. 
For each masing line, the ratio of the maximum opacity to t 22 is given, making use of the approximate treatment described in Sec.[3](numbers in 
pa rentheses) or with t he current method. The transitions in bold correspond to the maser transitions predicted by the present work and not predicted 
bv I Yates et al.1 dl997l) . The transitions with opacity ratio below 10~ 5 as well as the transitions that do not show population inversion are reported 
with the N symbol. 



Table 2. Same as Table[j]but forp-H 2 



Transition 


Frequency 


He 


P-H 2 


P-H 2 


o-H 2 


o-H 2 






(GHz) 


Green 


SET1 


SET2 


SET1 


SET2 


fad 






8.2(10.7) 


5.3 (7.1) 


9.1 (15.1) 


4.9 (6.8) 


7.6(10.1) 


3l,3 


- 2 2j o 


183 


0.64(1.22) 


0.98 (1.55) 


0.59 (1.29) 


1.08 (1.63) 


0.70(1.16) 


5 1.5 


-4 2 , 2 


325 


0.61 (1.24) 


0.85 (1.29) 


0.57 (1.35) 


0.93 (1.32) 


0.68 (1.22) 


7 5 ,3 


- 6 6 ,o 


437 


< 0.01 (0.02) 


0.16(0.23) 


0.09 (0.12) 


0.15 (0.22) 


0.13 (0.21) 


64.2 


-5 5 .i 


471 


0.46 (0.65) 


0.80(1.15) 


0.46 (0.70) 


0.83 (1.06) 


0.59 (0.97) 


53,3 


-4 4 ,o 


475 


0.49(1.40) 


0.77 (1.79) 


0.44 (1.42) 


0.86 (1.80) 


0.60(1.60) 


6 2 .4 


-7i, 7 


488 


N(N) 


<0.01 


NfN) 


<0.01 


<0.01 


2,., 


_ 2 , 2 


752 


N(N) 


N(N) 


N(N) 


N(N) 


N(N) 


9 2>8 


-83,5 


906 


0.45 (1.07) 


0.73 (1.41) 


0.43(1.00) 


0.79 (1.32) 


0.49 (0.98) 


4 2 .2 




916 


0.25 (0.56) 


0.36 (0.84) 


0.25 (0.73) 


0.39 (0.92) 


0.29 (0.81) 


5 2 .4 


-43,1 


970 


0.26 (0.80) 


0.55 (2.77) 


0.21 (0.81) 


0.60 (3.06) 


0.32(1.21) 


ll.l 


-0 ,o 


1113 


N(N) 


N(N) 


N(N) 


N(N) 


N(N) 


74.4 


-6 s ,i 


1173 


0.19 (0.41) 


0.41 (0.96) 


0.18(0.29) 


0.45 (0.98) 


0.23 (0.67) 


85,3 


-7 M 


1191 


0.18(0.21) 


0.15 (0.24) 


0.14 (0.18) 


0.19(0.29) 


0.24 (0.38) 


74.4 


- 8ij 


1345 


N(N) 


N(N) 


N(N) 


N(N) 


N(N) 


9 4 ,6 


- 10,., 


1435 


N(N) 


N(N) 


N(N) 


N(N) 


N(N) 


72.6 


- 63,3 


1441 


0.25 (0.41) 


0.49(1.61) 


0.23 (0.35) 


0.55 (1.88) 


0.29 (0.47) 


63.3 


_ 54,2 


1542 


0.15 (0.15) 


0.28 (0.35) 


0.13(0.11) 


0.32 (0.49) 


0.20 (0.35) 


73.5 


- 64, 2 


1766 


0.06 (0.07) 


0.24 (0.39) 


0.04 (0.06) 


0.30 (0.50) 


0.07 (0.11) 


84,4 


-7s,3 


2162 


< 0.01 (N) 


<0.01 


<0.01 


<0.01 


< 0.01 


93.7 


- 84,4 


2532 


<0.01 


0.16(0.18) 


< 0.01 (0.01) 


0.20 (0.23) 


< 0.01 


9 4 ,6 


-85,3 


2547 


N(N) 


<0.01 


<0.01 


<0.01 


N(N) 


6 2 .4 


-5 3 , 3 


2962 


N(N) 


< 0.01 


N(N) 


<0.01 


N(N) 


83,5 


~ 74_4 


3670 


N(N) 


< 0.01 


N(N) 


<0.01 


N(N) 



Notes. The ratios are still given as a function of the maximum 22 GHz opacity. 
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where v, is the rest frequency of line i, A(i) is the Einstein 
coefficient for spontaneous emission, g(i) is the degeneracy of 
the upper level of the transition and dra(z) is the population dif- 
fe rence per sub-leve l. This latter quantity is the one considered 
bv lYates et al.l 0997) and is defined as : 



d»(0 = P-^ 

\gl gu 



(21) 



where the indexes u and I stand respectively for the upper and 
lower level of transition i. In Eq. |20l the gas temperature Tj and 
turbulence velocity cr, are indexed according to the line, since 
in the comparison, we only retain the physical parameters that 
lead to the maximum opacity for this particular line. The opacity 
ratio thus correlat e to the population density ratio considered by 
lYates et al.l (1 19971) through : 



Tl 
T2 



<v 2 Y A(l) g(l)dn(l) 
v vJ A(2) g(2) dn(2) 



(22) 



2k B T l 2 




2k B T 2 2 


— +<n 




+ erf 


m 




m 



Since the temperatures and turbule nce of the models for 
which the ratios are given in Table 2 of lYates et all (1 1997b are 
unknown, we assume T\ - r 2 and cr\ = <x 2 while performing the 
conversion. The maximum error introduced by this assumption 
is typically of a factor ~ 2-3 (considering that the 22 GHz has 
its maximum inversion at 600 K as stated in Yates (1997). ) By 
considering col. 5 of Table |4] it can be s een that the value s cur- 
rently obtained and the one obtained by lYates et all (1 19971) can 
differ by factors larger than 2-3. The d ifferences m i ght arise from 
the use of plane-parallel geometry in lYates et all (1 19971) against 
spherical geometry in the current work, or because we do not 
consider the effect introduced by the pumping by dust photons 
in the current study. By examining the results of the comparison 
reported for 0-H2O in Tableland for the lines with frequencies 
higher than 1.3 THz (A < 230 yum), we can see that the ratios 
can take high values above this threshold frequency and tend to 
increase while the frequency increases. Since this behaviour is 
correlated with the shape of the dust emission, it can be guessed 
that a quantitative comparison between the two works is plagued 
by the absence of consideration of dust radiation in the present 
work. 

By compari s on to prev ious calculations 

dNeufeld & Melnickl I199U lYates et all 119971) . we find that 
the 85,4 - 7 6 ,i (1165 GHz), 7 4 , 4 - 65,1 (1173 GHz), 8 5;3 - 7 6>2 
(1191 GHz) lines can show population inversion with substantial 
opacities, independently of the set of rate coefficients used. 
These lines were not predicted to be masers in the previous 
studies. The rest of th e lines con s idered in this work agree 
with the predictions of lYates et al.l d 19971) with respect to the 
possibility of masing action, at least for the lines that show 
substantial inversion. 

Considering the effects introduced by the collisional rate co- 
efficients, we find that the variations from one set to another are 
moderate. The lines which are found to be inverted are the same, 
independently of the set considered, and most of the opacities 
show variations of the order of a factor 2 or less. Our method 
to treat the RT of the maser lines introduces large differences in 
the computed line opacities; indeed, for most of the lines, the 
opacities derived with the exact treatment are lower by up to a 
factor 4-5 in comparison to the opacities derived with the ap- 
proximated treatment. 



Interestingly, we find that many lines can show substantial 
inversion. Th i s is in agreement with the observations made by 
iMenten et al.l (|2008) who report the fluxes of 9 maser transitions 
of water. Note that the 437 GHz 75 3 - 65,0 whose first detection 
is reported in the latter study is expected to show a substantial 
inversion (as was found in Yates et al., 1997) wh ile thi s line was 
not found to be inverted by Neufeld & Melnick| (Il99ll) . 



Finally, Humphreys et al.l d2001l) presented models for the 
excitation of the 22, 183, 321 and 325 GHz water lines in the cir- 
cumstellar envelopes (CSE) of AGB stars. They found the differ- 
ent masing lines to be excited in specific regions of the CSE, thus 
predicting different morphologies of the emission maps for the 
masing lines they considered. A direct comparison with their re- 
sults is not possible, since in their study, the model that describes 
the circumstellar envelope is rather elaborated. Indeed, they first 
compute the H 2 volume density, gas temperature and expansion 
velocity from an hydrodynamical code. As a second step, they 
randomly define masing spots in the envelope and solve the ID 
transfer problem for the various spots using the LVG approxi - 
mation. The physical parameters that concern each emitting spot 
is defined according to its distance to the star with respect to the 
parameters obtained from the hydrodynamical model. Finally, 
the level populations derived from the LVG calculation are then 
used in a code that propagates the radiation and that takes into 
account maser saturation in order to calculate the emission of 
each spot. Performing such a model is out of the scope of the 
current study, but it is possible to qualitatively discuss their re- 
sults with the current ones. In Fig. [9] we show the opacities of 
various water maser lines as a function of the H 2 volume density. 
The rate coefficients used consider o-H 2 as a collisional partner 
and are taken from SET2. The opacities are obtained using the 
treatment for maser propagation and for each line, the opacity is 
given at 500K and 1000K. 



One of the conclusion drawn by Humph reys et al.l (12001 1) is 
that the 321 and 22 GHz are sensitive to the acceleration zone 
of the CSE. Additionally, the 321 GHz line is predicted to trace 
only a region close to the star while the 22 GHz emission ex - 
tends further outside (see Fig. 9 of [Humphre ys et al.l 12001 1) . 
Considering the opacities reported in Fig. [9] we can see that 
the current results support those findings. Indeed, at densities 
~ 10 9 cm -3 , the 22 GHz line is the only masing line which is 
still subsequently inverted. At such a high H 2 density, the other 
masing lines are found to be collisionaly quenched. Moreover, 
the 321 GHz line is inverted only at relatively high H 2 den- 
sities, ie. n(H 2 ) > 2 10 7 cm~ 3 which support the fact that this 
line on ly traces the innermost part of the CSE. Humphreys et al. 
(2001) also find that the emission in the 325 GHz line resem- 
bles that of the 22 GHz. In the current study, by examining the 
behaviour of the opacities of those two lines shown in Fig. [9] 
we can indeed see t hat these lines have sim ilar excitation con- 
ditions. Additionaly Humphreys et al. (2001) predicted the 183 
GHz line to be inverted in the outermost part of the CSE. In this 
region, the H 2 density is too low to produce inversion in the other 
masers. Referring to Fig.|9] we can see that at n(H 2 ) ~ 10 6 cm" 3 , 
the 183 GHz line is indeed subsequently inverted while the other 
lines are not. Finally, considering the behaviour of the opacities 
reported in Fig. [9] it is possible to qual itatively predict the mor- 
phology that the lines not considered by Humphr eys et al.l (1200 ll) 
should have. As an example, the emission in the 380 GHz line 
should resemble the 183 GHz one and the 439 and 470 GHz lines 
should behave similarly to the 321 GHz line. 
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Table 3. Comparison of the current results with lYates et al.l(ll997l) 's results. 



o/p-H 2 Transition 


Frequency (GHz) 


Yates97 


current 


Yates97/current 


6l,6 _ 52.3 


22 


1.0 


1.00 


1.00 


102,9 - 9 3 .6 


321 


0.59 


0.88 


0.67 


4,, 4 - 3 2a 


380 


1.46 


1.78 


0.82 


64.3 - 55.o 


439 


0.13 


0.63 


0.20 


75.2 - 6 6 ,i 


443 


0.14 


0.05 


2.89 


4 2 ,3 - 33,0 


448 


1.12 


1.49 


0.75 


5 3 ,2 - 4 4 ,i 


621 


0.36 


1.69 


0.21 


63,4 - 54,1 


1158 


0.89 


0.80 


1.11 


74.3 _ 65.2 


1278 


0.09 


0.49 


0.19 


82.7 _ 73.4 


1296 


1.74 


0.76 


2.29 


62.5 _ 53.2 


1322 


3.11 


0.34 


9.17 


84,5 _ 7 52 


1885 


0.31 


0.02 


15.7 


5 2 ,3 - 4 3 , 2 


1918 


1.74 


< 0.01 


> 173 


83.6 - 74.3 


2245 


0.86 


< 0.01 


> 85 


73.4 - 64.3 


2567 


1.17 


< 0.01 


> 116 


3 1.3 - 22,0 


183 


0.83 


1.22 


0.68 


5l,5 -4 2 . 2 


325 


1.02 


1.24 


0.82 


75,3 - 66.0 


437 


0.10 


0.02 


4.80 


64,2 _ *5,i 


471 


0.07 


0.65 


0.11 


53,3 - 44,0 


475 


0.25 


1.40 


0.18 


92,8 - 83,5 


906 


0.60 


1.07 


0.55 


42,2 - 33,1 


916 


0.48 


0.56 


0.86 


5 2 ,4 - 4 3 ,i 


970 


1.19 


0.80 


1.49 


72.6 - 63,3 


1441 


1.07 


0.41 


2.62 


63,3 _ 54,2 


1542 


0.11 


0.15 


0.72 


73.5 - 64.2 


1766 


0.15 


0.07 


2.21 



Notes. Comparison of the opacities obtained in the present work (column 4) to the one estimated from Yates et al. ( 1997) (column 3). In this latter 
case, the opacities are obtained from the population density ratio they reported and making use of Eq. f22](see text for details). The last column 
gives the opacity ratio obtained in the two studies. 



5. Conclusions 

The present work deals with obtaining a self-consistent solution 
for the problem of maser propagation within the scope of non- 
local radiative transfer. Test cases are presented for the water 
molecule in which the importance of dealing with maser prop- 
agation is emphasised, by comparison to a simplified treatment 
for the radiative transfer. The main conclusions are : 

1 . A self-consistent solution for the radiative transfer in masing 
lines can be obtained within the scope of the short character- 
istic method, with subsequent modification of the algorithm 
used to treat thermal lines due to the integro-differential na- 
ture of the transfer equation. 

2. The preconditioning of the statistical equilibrium equations 
is not feasible using the current method and the Ng acceler- 
ation technique cannot be used, so that the algorithm has the 
convergence rate of the Lambda iterative scheme. 

3. The test cases performed for the water molecule have shown 
that neglecting maser propagation can lead to substantial er- 
rors in the estimate of the masing line opacities, with differ- 
ences greater than a factor 2. Dealing with strongly inverted 
masers, it means errors of several orders of magnitude in the 
prediction of line intensities. 

4. In the region where strong masing lines are present, the ra- 
diation field in these lines perturbs the whole population of 
the molecule. Therefore, including maser propagation also 
modifies the determination of line opacities for purely ther- 
mal lines and in the present test cases, maximum differences 
of the order of a few 10% are found between the two meth- 
ods adopted. Nevertheless, the number of lines affected is 
rather low and most of the lines have identical optical depths 



in the two treatments. Given the larger computational time 
required to solve the maser propagation, the approximation 
used in this study is a good alternative if one is concerned 
with interpreting the emission of non-masing lines. 
5 . We give predictions of H2O maser opacities for physical con- 
ditions typical of circumstellar envelopes of AGB stars. In 
addition, the predictions based on various sets of rate coef- 
ficients involving He, p-I-b or 0-H2 are compared. We find 
that from one set to another, the relative opacities of the mas- 
ing lines are similar. However, the absolute scale of the opac- 
ities are within a factor ~ 2 depending on the set used. 
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Fig. 9. Opacity as a function of H2 volume density, for various 0-H2O and P-H2O lines which are found to be substantially inverted 
at T K = 500 K (red lines) and TV = 1000 K (black lines). 
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